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SUMMARY 


The formation of nitrogen oxides (N0 X ) during combustion of methane, pro- 
pane, and a jet fuel, JP-4, was investigated in a jet-stirred combustor. The 
results of the experiments were interpreted using reaction models in which the 
nitric oxide (NO) forming reactions were coupled to the appropriate hydrocarbon 
combustion reaction mechanisms. Comparison between the experimental data and 
the model predictions reveals that the CH + N 2 reaction process has a signif- 
icant effect on nitric oxide formation especially in stoichiometric and fuel- 
rich mixtures. Reaction models were assembled that predicted nitric oxide 
levels that were in reasonable agreement with the jet-stirred combustor data 
and with data obtained from a high-pressure (5.9 atm (0.6 MPa)), prevaporized, 
premixed, flame- tube- type combustor. The results of the experiments and the 
theoretical studies also suggested that the behavior of hydrocarbon mixtures, 
like JP-4, may not be significantly different from that of pure hydrocarbons. 
Application of the propane combustion and nitric oxide formation model to the 
analysis of NO x emission data reported for various aircraft gas turbines showed 
the contribution of the various nitric oxide forming processes to the total 
NO x formed. 


INTRODUCTION 

The ability to describe the kinetics of nitrogen oxide formation during 
combustion of hydrocarbon-based fuels is an important requirement for the 
development of analytical combustor models such as those for gas turbines. The 
mechanism of nitric oxide formation from atmospheric nitrogen has been studied 
extensively, and during combustion of fuel-lean and near-stoichiometric fuel- 
air mixtures, the formation of nitric oxide is generally understood by what is 
often called the extended Zeldovich reaction mechanism: 


O + N 2 NO + N 
N + 0 2 -* NO + O 
N + OH + NO + H 


This mechanism, when coupled with reactions that describe the oxidation of hydro- 
carbon species and the formation of atomic oxygen and the hydroxyl radical, pre- 
dicts nitric oxide production rates and concentrations that are a function of 
both temperature and equivalence ratio. However, these predictions have not 
always been in agreement with observed nitric oxide levels, particularly for 
hydrocar bon- air flames. For example, during an investigation of the formation 
of nitrogen oxides in hydrogen, carbon monoxide, and propane flames within a 
jet-stirred combustor, Engleman et al. (ref. 1) compared experimental nitric 
oxide concentrations with those predicted from a chemical kinetic model of the 



appropriate system. The extended Zeldovich mechanism adequately modeled forma- 
tion of nitric oxide in the nonhydrocarbon (H 2 and CO) systems at all equiva- 
lence ratios. However, in the propane flames, the observed levels were signif- 
icantly underestimated. More recently, Engleman et al. (ref. 2) and Wakelyn 
et al. (ref. 3) reported similar differences between observed and theoretical 
nitric oxide levels when the extended Zeldovich reaction scheme is the only path 
for nitric oxide formation. 

Explanations have been proposed to account for these differences between 
observed and predicted nitric oxide levels in hydrocarbon-air flames. It has 
been suggested (refs. 4 to 6) that excess amounts ("super-equilibrium") of 
atomic oxygen and the hydroxyl radical are responsible for the higher observed 
nitric oxide levels and that the current hydrocarbon oxidation mechanisms do 
not adequately predict these 0 and OH concentrations. Malte et al. (ref. 7) , 
however, measured hydroxyl radical and atomic oxygen concentrations formed dur- 
ing combustion of methane in a jet-stirred combustor and found that these con- 
centrations were not high enough to account for the observed NO concentrations. 
They concluded that the lack of agreement between experiment and theory indi- 
cates a possible influence of reactions between hydrocarbon fragments and mole- 
cular nitrogen leading to nitric oxide, as first suggested by Fenimore (ref. 8) . 

Fenimore proposed that the reaction between the CH radical and molecular 
nitrogen 


CH + N 2 HCN + N 


produced intermediates that were subsequently oxidized to nitric oxide. Subse- 
quent experimental studies (refs. 9 to 1 1 ) have verified the presence of HCN 
in various hydrocarbon-air flames. Thus, the oxidation of HCN and N could be 
a significant source of NO in flames. Engleman et al. (ref. 2) and Wakelyn et 
al. (ref. 3) have had some success in predicting the NO levels observed during 
combustion of methane and propane in a jet-stirred combustor. The agreement 
between observed and calculated NO levels was reached by adjusting the rate coef- 
ficient for the reaction CH + N 2 * HCN + N. The resulting rate coefficients 
were different for the two studies and furthermore were not in good agreement 
with the recent rate coefficient expressions reported by Matsui and Nomaguchi 
(ref. 12) and Blauwens et al. (ref. 13). 

The objective of the study reported here was to experimentally and theoret- 
ically investigate the production of nitrogen oxides during combustion of vari- 
ous hydrocarbon fuels, including a jet fuel, and attempt to assemble a chemical 
kinetic combustion and nitric oxide formation model that is consistent with 
experimental results and recent kinetic data. The experiments were performed 
in a jet-stirred combustor and the fuels studied were methane, propane, and the 
jet fuel, JP-4, at fuel-air equivalence ratios from 0.7 to 1.4. Since propane 
is often used as a surrogate jet fuel in many combustion experiments, a propane- 
air experiment was carried out under the same test conditions as the jet-fuel- 
air experiments to provide a direct comparison of the NO levels at identical 
mass loadings. Such a comparison should indicate whether a propane combustion 
and NO formation model would also be a good surrogate model for a jet fuel. 
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EXPERIMENTAL APPARATUS AND PROCEDURE 


The combustion device used in these studies consisted of a jet-stirred 
combustor identical to the device described by Wakelyn et al. (ref. 3) with the 
exception that the injector , as well as the shell, was made of castable zirconia 
(see fig. 1). Mixtures of fuel and air were fed to a hemispherical injector 
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Figure 1.- Cast zirconia jet-stirred combustor. 
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containing 40 evenly distributed holes of 0.051 -cm diameter that opened into 
a hemispherical combustion cavity of 12.7-cm^ volume. Combustion products were 
exhausted through 25 holes of 0.318-cm diameter that were evenly distributed 
over the hemispherical cap of the combustor. Two additional holes of 0.41 6-cm 
diameter were provided in the cap for use as gas sampling and thermocouple 
ports; these holes were positioned on opposite sides of the cap along a horizon- 
tal line located sufficiently above the injector surface to permit probes to 
pass over the injector without interference. 

The change to an all-zirconia combustor design served to avert flashback 
during combustion of liquid hydrocarbon fuels. The former design, with an 
inconel injector, allowed sufficient heat to be transferred through the injector 
wall to raise the temperature of the air-fuel mixture within the injector to 
its autoignition point, resulting in combustion within the injector itself. 

This condition could not be avoided for fuel-air equivalence ratios from 0.9 
to 1.2 without resorting to intentional cooling of the combustion zone (e.g., 
Singh (ref. 14) employed a nitrogen diluent). Zirconia provided sufficient 
resistance to heat transfer across the injector wall to prevent flashback. 

The combustion of liquid fuels in this device required a different approach 
for premixing and metering fuel and air than that described in reference 3 for 
propane. Since experimental limitations did not readily allow the use of heat 
alone to effect fuel evaporation, the approach taken by Wright (ref. 15) was 
followed in which the fuel was both heated and atomized. In the present work, 
the fuel was finely atomized at approximately 190°C. It was expected that this 
technique would suffice for JP-4 provided that atomization occurred sufficiently 
near the combustor injector to prevent condensation of the high-boiling-point 
fractions of the fuel at temperature below the 100-percent evaporation point. 

This expectation appeared to be supported by the calculations of reference 16, 
which indicate that no residue from an average refinery run of JP-4 would 
be lost through condensation. Results of subsequent experiments with this test 
configuration confirmed these calculations, since no high-boiling-point fractions 
could be observed inside the test apparatus after more than 20 liters of JP-4 
had passed through the system during the combustion experiments. 

The apparatus used to vaporize the fuel and control the flow of fuel-air 
mixtures into the combustor is depicted in figure 2. Fuel reservoirs consisted 
of closed vessels pressurized with nitrogen gas. Fuel flow rate was measured 
by rotameters and controlled by fine-metering valves. Valves located downstream 
of the rotameters kept sufficient pressure on the fuel to avoid bubble formation 
from dissolved nitrogen as it passed through the rotameters. Mixtures of dif- 
ferent fuels could be obtained by adjusting flow rates from individual fuel 
reservoirs. 

It was found that changes in room temperature could change fuel density 
and viscosity sufficiently to cause flow rate errors of as much as 1 percent 
per degree Celsius. Since room temperature ranged from 18°C to 27°C, a con- 
trolled temperature water bath was employed to keep fuel passing through the 
rotameters at a temperature of 40°C ± 0.5°C. 
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Figure 2.- Schematic diagram of combustor fuel and air supply system. 


An internal-mix, two-fluid nozzle was used for atomization in which the 
flow rate of impinging nitrogen gas was maintained at a constant 20 L/min. This 
provided fine atomization over the entire range of fuel flows, 0.5 mL/min to 
7.5 mL/min. Care was taken in the nozzle design, and in the orientation of the 
expansion chamber into which the spray was directed, to prevent fuel from drip- 
ping because this invariably resulted in unstable combustion. 

Downstream of the expansion chamber, oxygen and nitrogen were mixed with 
the fuel- vapor-nitrogen mixture exiting from the expansion chamber. Flow rates 
of oxygen and nitrogen were controlled by fine-metering valves and monitored by 
mass flowmeters. Total oxygen flow rates ranged to 13 L/min and total nitrogen 
flow rates ranged to 50 L/min. Preheat temperature of the entire combustible 
fuel-air mixture was measured within the combustor stem with a platinum/ 
platinum+1 4-percent-rhodium thermocouple at a point sufficiently distant from 
the combustion cavity to avoid temperature feedback effects from the combustion. 
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Combustion cavity temperature was measured and gas samples were taken in 
accordance with the procedures described in reference 3. A traversing iridium/ 
indium+40-percent-rhodium thermocouple measured flame temperature; and a water- 
cooled, aerodynamically quenched, sampling probe of stainless steel withdrew 
gas samples for analyses of NO and N0 X with a chemiluminescence analyzer. Gas 
samples were taken at a probe depth into the combustion cavity of 0.15 to 
0.25 cm. The only difference between the sampling procedures in the present 
work and those described in reference 3 was the sampling dilution ratio: 
whereas a dilution ratio of 6 to 1 was satisfactory for propane and JP-4, a 
ratio of approximately 8 to 1 was required for methane combustion to prevent 
water condensation in the sampling lines. 

Experiments were performed at ambient (atmospheric) pressure at a combustor 
inlet temperature of 27°C (300 K) with the methane-air mixtures, 27°C and 182°C 
(455 K) with the propane-air mixtures, and 182°C with the JP-4-air mixtures. 

The combustor mass loadings were between 0.072 and 0.074 g/cm^-sec correspond- 
ing to residence times between 2.4 and 2.6 msec. The fuel-air equivalence ratio 
was varied between 0.7 and 1.4. 


CHEMICAL KINETIC MODELS 

The reaction models that were assembled for comparison with the experimen- 
tal data are, for the most part, based on chemical kinetic and mechanistic infor- 
mation found in the literature. It was not the intent, however, to assemble an 
extensive list of reactions, but rather to develop reaction mechanisms contain- 
ing enough detail to describe essential features of the oxidation and nitric 
oxide formation. The rate coefficients selected for the reactions in the vari- 
ous mechanisms were obtained from the literature whenever possible. (See 
refs. 3, 10, 12, and 17 to 37.) It must be emphasized at this point that the 
proposed mechanisms are not necessarily complete since very little kinetic and 
mechanistic information is available on the hydrocarbon-oxygen systems, espe- 
cially on reactions involving the lower molecular weight components such as CH 
and CH2 that are believed to be important species in the production of nitric 
oxide in fuel-rich mixtures (ref. 8) . 

The methane oxidation mechanism that was assembled for this study is listed 
in table 1 . The high- temperature combustion of methane has been investigated 
for many years, and most of the basic features of the combustion mechanism are 
known. The methane combustion mechanism used in this study contains, we believe, 
all these essential features. 

The propane combustion mechanism used in this study is listed in table 2. 
This mechanism is based on the results of an experimental and theoretical study 
of propane oxidation behind shock waves (ref. 30) . The propane combustion mech- 
anism is a rather simplified mechanism in which propane is assumed to rapidly 
decompose at combustion temperatures and react to form lower molecular weight 
hydrocarbons such as the methyl radical, ethylene, and acetylene, which are sub- 
sequently oxidized. A minimal number of reactions were used to describe the 
transformation of propane to the lower molecular weight fragments. The mech- 
anism presented in table 2 distinguishes between the normal propyl and isopropyl 
radicals. This distinction was required to explain the shock tube experiments. 
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but was unnecessary for the present study. By using an "average" rate coeffi- 
cient for reactions involving the propyl radical, the propane mechanism can be 
simplified even further without affecting the calculated nitric oxide levels. 

The complete propane combustion scheme used in this study consisted of the 
reactions in table 2 combined with the reactions in table 1 . 

As noted previously, Fenimore (ref. 8) suggested that the reaction between 
the CH radical and molecular nitrogen could produce the intermediates HCN and 
N that eventually were transformed to nitric oxide. The reactions that were 
selected to describe the transformation of HCN and N to nitric oxide, as well 
as the formation and consumption of the CH radical, are listed in table 3. 

Also included in table 3 are reactions describing the formation of NO through 
N 2 O which contributes to NO formation during the combustion of very fuel-lean 
mixtures. The importance of these reactions in determining nitric oxide levels 
is discussed later. The reactions listed in table 3 combined with the methane 
and propane oxidation schemes listed in tables 1 and 2 formed the kinetic models 
used in this study to calculate nitric oxide levels. 

To model the experiments, a computer program was used which describes the 
jet-stirred combustor as a perfectly stirred reactor (PSR) in steady-state oper- 
ation. The program is basically a combination of the equilibrium program of 
Gordon and McBride (ref. 38) and the PSR solution algorithm of Jones and 
Prothero (ref. 39) . Inputs to the program include a reaction mechanism, initial 
composition, temperature and pressure, and a heat loss parameter to account for 
the nonadiabatic operation of the reactor. Solutions are provided at desired 
mass loadings. 


RESULTS AND DISCUSSION 

Comparison Between Experimental Results and Chemical Kinetic Models 

The results of the methane-air, propane-air, and the JP-4-air experiments 
are plotted in figures 3 to 6 where the measured reaction temperature and nitro- 
gen oxide (N0 X ) concentration in parts per million (ppm) are plotted against 
the fuel-air equivalence ratio. The data were obtained at mass loadings between 
0.072 and 0.074 g/cm^-sec which correspond to residence times between 2.4 and 
2.6 msec. Maximum nitrogen oxide concentrations occurred in the range of equi- 
valence ratios between 1.0 to 1.1. The NO x concentrations and flame temperature 
for the propane-air and JP-4-air experiments (see figs. 5 and 6) were essen- 
tially the same for the stoichiometric and fuel-lean mixtures. For the fuel- 
rich mixtures, the observed NO x levels were higher for the jet-fuel experiments; 
however, as shown later, this can be attributed to higher flame temperatures. 

The ability of the proposed mechanisms to predict the experimentally mea- 
sured nitrogen oxide levels was assessed by simulating the jet-stirred combustor 
experiments using the PSR program and the reactions and rate coefficients listed 
in tables 1 to 3. These simulations were made for fuel-air equivalence ratios 
between 0.8 and 1.3. The heat loss parameter in the PSR program was adjusted 
to provide a reasonable match between the calculated and measured temperatures 
for a mass loading of 0.073 g/cm 3 -sec. The flame temperatures used in these 
simulations are noted in figures 3 to 6. 
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Figure 3.- Experimental and calculated results for methane-air mixtures 
with combustor inlet temperature of 300 K. Calculations made with 
k 66 = 8 x 10 10 exp (-6844/T) . 
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Figure 5.- Experimental and calculated results for propane-air mixtures 
with combustor inlet temperature of 455 K. Calculations made with 
k 66 = 8 x 10 10 exp (-6844/T) . 
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Figure 6.- Experimental and calculated results for JP-4-air mixtures with 
combustor inlet temperature of 455 K. Calculations made with propane 
model with = 8 x 10 10 exp (-6844/T) • 
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The results of these simulations revealed that the assembled chemical 
kinetic models overpredicted the measured nitric oxide levels by a factor of 
up to 4. To better understand the proposed combustion and nitric oxide forma- 
tion models in relation to the experimental data, a parametric study was per- 
formed in which the rate coefficients for selected reactions were varied and 
the effect on the nitric oxide levels was noted. Such a study would hopefully 
identify key reactions which govern the formation of NO and determine whether 
it was possible to achieve good agreement between experiments and computations 
by making reasonable adjustments to the model. Selected results from the study 
are given in table 4 for simulations of the methane experiment and one of the 
propane experiments (at temperature T of 455 K) . 

The results of the parametric studies revealed that the reactions which 
describe the formation of NO through the oxidation of HCN and N are needed in 
the mechanism to account for the observed N0 X concentrations. When the only 
path for nitric oxide formation is through the extended Zeldovich reactions, 
the model underpredicts the observed N0 X levels especially for the stiochio- 
metric and fuel-rich mixtures (see table 4) . 

The reaction between the CH radical and molecular nitrogen (reaction 66 ) 
has the greatest effect on the calculated nitric oxide concentrations, as 
expected. In fact, by adjusting the rate coefficient for this reaction to 
k 66 = 8 x TO 10 exp(-6844/T) , the nitric oxide levels predicted by the kinetic 
models are in excellent agreement with the experimental results. The adjusted 
rate coefficient expression yields rate coefficients that are considerably 
smaller than the values obtained from the expressions reported in references 12 
and 13. The rate coefficient reported in reference 12, kg£ = 4 x 10 11 
exp(-6844/T) , is a rough estimate based on experimental data obtained at 2500 K, 
and the uncertainty associated with this estimate was not reported. The rate 
coefficient reported in reference 13, kg^ = 8 x lO ^ 1 exp (-5536/T) , represents, 
according to the investigators, an upper limit with an uncertainty of at least 
a factor of 2. On the basis of this information, it is difficult to determine 
whether the adjusted rate coefficient is in reasonable agreement with these 
expressions. However, the adjusted rate coefficient expression is reasonable 
for a biomolecular reaction of the type represented by reaction ( 66 ) . Further- 
more, the expression is in fair agreement with a theoretically calculated 
expression (ref. 2), k 55 = 1 x 10 11 exp(-9562/T) , and the expression, 
k 66 85 1-5 x 10 11 exp (-9562/T) , reported in reference 3. 

As the results presented in table 4 indicate, the other reactions that 
had a large influence on the calculated NO levels were those that controlled 
the formation and destruction of the CH radical. In the methane combustion model, 
the major sources of the CH radical are reactions (59), (60), and (61) with some 
contribution from reaction (23) for fuel-rich mixtures (equivalence ratios 
4> of approximately 1.3): 


(59) 

CH 2 + H 

■+ CH + H 2 

(60) 

ch 2 + 0 

■* CH + OH 
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(61) CH 2 + OH > CH + H 2 0 
(23) C 2 H +0 + CH + CO 

The major paths for CH destruction are reactions (25) and (62) with some influ- 
ence of reactions (63) and (64) for very rich mixtures (4> m 1.3): 

(25) CH + 0 2 -► HCO + 0 

(62) CH + C0 2 -*■ HCO + CO 

(63) CH + H -► C + H 2 

(64) CH + OH * C + H 2 0 


Reaction (25) has a large effect for fuel-lean and near-stoichiometric mixtures, 
while reaction (62) has a large effect over the entire range of equivalence 
ratios examined. 

In the propane combustion model, the major sources of the CH radical are 
again reactions (59), (60), and (61), with a much greater influence by reac- 
tion (23) than in the methane oxidation scheme. The increased effect of reac- 
tion (23) occurs because a primary path in the propane oxidation mechanism is 
the formation of acetylene which is the source of the C 2 H radical, while in the 
methane oxidation mechanism, oxidation of the methyl radical is the primary path. 
Acetylene is produced in the methane mechanism only after the formation of C 2 Hg. 
The major paths for CH destruction in the propane model are reactions (25) 
and (62) . 

As note previously, an objective of the parametric study was to determine 
whether good agreement between experimental results and the kinetic models 
could be achieved by making "reasonable" adjustments to the rate coefficients 
for key reactions. The adjustments were focused mainly on the reactions which 
controlled the formation and destruction of the CH radical. Since the assembled 
chemical kinetic models over-predicted the observed NO x concentrations, a 
decrease in the nitric oxide levels could be produced by decreasing the rate 
of production of the CH radical and/or increasing its rate of consumption. As 
noted previously, the primary sources of the CH radical are different for the 
methane and propane models. In the methane mechanism, the rate of reaction (61) 
has the largest influence on CH production, while in the propane model, the rate 
of reactions (61) and (23) control the formation of CH. Attempts to decrease 
the rate coefficients of these reactions to achieve reasonable agreement between 
the calculated and experimental results were not very successful. While reason- 
able agreement could be achieved separately for either the methane or the pro- 
pane results, a satisfactory fit to both simultaneously could not be obtained. 

Variations in the rate coefficients of the CH consuming reactions revealed 
that the rate of reaction (62) had a significant effect on calculated NO levels 
at all equivalence ratios, while reactions (25) , (63) , and (64) had largest 


13 



effects for stoichiometric and fuel-lean mixtures. Attempts to achieve agree- 
ment between calculated and experimental results by increasing the rate coeffi- 
cients for reactions (25), (63), and (64), either separately or together, did 
not prove successful because of the dominant effect of these reactions in 
the lean and stoichiometric mixtures. For example, when good agreement was 
obtained in fuel-rich mixtures, the models underpredicted the observed NO 
levels in the lean and stoichiometric mixtures. Increasing the rate coefficient 
of reaction (62), however, was more successful. Very good agreement was 
achieved by adjusting the rate coefficient of reaction (62) to 
k 62 = ** x 10 12 T 1 / 2 exp(-3020/T) which is a factor of 100 larger than the 
initial expression. Such a large adjustment may not be reasonable, however, 
since the resulting rate coefficient would have a preexponential factor 
(1 X 10 12 T 1 / 2 ) that is unusually large for a reaction between molecular 
species. 

In summary, the results of the parametric study revealed that the best 
agreement between the predictions of the chemical kinetic models and the 
stirred reactor data for methane and propane was achieved by adjusting the 
rate coefficient for reaction (62) or reaction (66) . Since the predicted 
nitric oxide levels were most sensitive to the rate of reaction (66) and 
since the required adjustment to the rate coefficient for reaction (62) was 
considerable and somewhat unreasonable, it was concluded that the most reason- 
able and kinetically consistent model was achieved when the rate coefficient 
for reaction (66) was set to 8 x 10^ exp (-6844/T) . The nitric oxide levels 
(expressed as ppm NO x ) predicted by the methane and propane models using the 
adjusted rate coefficient for reaction (66) are compared with the experimental 
results in figures 3 to 6. Note that the predictions from the propane model 
are also in good agreement with the JP-4 experimental results (fig. 6) . 
Apparently, the difference in the NO x levels that were observed between the 
propane and JP-4 experiments can be attributed to the difference in flame 
temperatures. The good agreement between the predictions of the propane model 
and the JP-4 experimental results also suggests that the propane combustion and 
NO formation model can be used as a surrogate kinetic model for a jet fuel. 


Comparison Between Chemical Kinetic Models and Other Data 

As an additional check of the propane combustion and NO formation model, 

N0 X concentrations predicted by the model were compared with N0 X emissions 
reported by Anderson (ref. 40) for premixed, prevaporized propane-air mixtures 
in the NASA Lewis Research Center 1 0-cm-diameter flame- tube- type combustor. 

These experiments were conducted at an inlet temperature of 590 K and a pressure 
of 5.9 atm (0.6 MPa) with a reference velocity of 23 m/sec. The concentration 
of NO x was measured at a point 46 cm downstream of the flameholder for a span 
of fuel-air equivalence ratios ranging from the lean blowout limit to slightly 
richer than stoichiometric. The results of these studies, given in figure 7, 
show the strong dependence on equivalence ratio that can be expected from 
premixed, prevaporized operation. The nitric oxide levels predicted by the 
propane model (with the adjusted rate coefficient for reaction (66)) are also 
presented in figure 7. These calculations were made assuming adiabatic opera- 
tion (ref. 40) . Since the N0 X measurements were made at the combustor center- 
line, the adiabatic assumption was not unreasonable. The agreement between the 
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analytical model and the data is quite good. Also shown in figure 7 are the 
contributions that the extended Zeldovich mechanism, the N 2 O reaction mechanism, 
and the CH + N 2 reaction scheme each make to the total N0 X . The N 2 O reactions 
contributed to the total N0 X levels at equivalence ratios less than 0.7, while 

100 r- 



Fuel -air equivalence ratio 

Figure 7.- Experimental NO x levels for propane-air mixtures in flame- tube- 
type combustor (ref. 40) compared with NO x levels predicted by propane 
combustion and NO formation model. 


the oxidation of the HCN and N produced by reaction (66) contributes to the 
total NO x when the equivalence ratio is greater than 0.8. 


Analysis of Aircraft NO x Emission Data 

Since the propane combustion and nitric oxide formation kinetic model pro- 
vided a reasonably good description of the JP-4 experimental results, the model 
was used to investigate the combustion kinetics and the nitric oxide formation 
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processes for typical aircraft gas turbine operating conditions. This was 
accomplished by comparing the nitric oxide levels predicted by the model with 
the N0 X emissions for various aircraft gas turbines as originally presented by 
Lipfert (ref. 41). Lipfert demonstrated that N0 X emission data, expressed in 
g NC> 2 /kg fuel and adjusted to a common humidity level (O.Olg ^O/g dry air), 
correlated well with combustor inlet temperature for a wide variety of engines. 
Figure 8 reproduces the Lipfert correlation. 



Figure 8.- Comparison of predictions from kinetic models with Lipfert data 
(ref. 41). In the complete mechanism, kgg = 8 x 10 10 exp (-6844/T) . 


To compare the model with the Lipfert data, it was necessary to determine 
a representative combustor residence time. The procedure reported by Blazowski 
et al. (ref. 42) was used to calculate the residence time. Blazowski et al. 
computed a combustor residence time by comparing the N0 X levels predicted by 
a model with the Lipfert data for a representative condition - a combustor inlet 
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temperature of 675 K, a humidity of O.Olg H 2 0/g dry air, and an emission index 
of 15g N0 2 /kg fuel. These conditions correspond to sea level static operation 
at 300 K, a pressure ratio of 13.8, and a representative compressor efficiency 
of 0.85. Having determined the residence time for this condition, it was 
assumed to be valid for all other inlet conditions in the Lipfert correlation 
for a static temperature of 300 K and a compressor efficiency of 0.85. 

With these representative operating conditions inserted into the PSR model 
and adiabatic operation with a stoichiometric fuel-air mixture assumed, the 
representative combustor residence time was computed to be approximately 
0.8 msec. By using this residence time in PSR calculations for the other inlet 
conditions in the Lipfert correlation, the solid curve shown in figure 8 was 
generated. 

Excellent agreement with the Lipfert data was obtained throughout the 
entire range of inlet conditions. Also shown in figure 8 are curves represent- 
ing the NO x levels predicted by the model when the CH + N 2 reaction scheme is 
removed and by a quasi-chemical kinetic model. In the quasi-chemical kinetic 
model, the fuel-oxygen system is assumed to be in chemical equilibrium at the 
adiabatic flame temperature and the formation of nitric oxide is controlled by 
the kinetics of the extended Zeldovich reaction scheme. A comparison of the 
three curves reveals the contribution of several factors to the total NO x . 

The difference between the curves generated by the complete model and the model 
with only the extended Zeldovich and N 2 0 reaction schemes represents the con- 
tribution of the CH + N 2 reaction scheme. This contribution decreases as the 
temperature increases. The influence of the N 2 0 reaction scheme on NO x levels 
was found to be negligible for the conditions in this analysis. 

The difference between the curves generated by the model containing only 
the extended Zeldovich and N 2 0 reaction schemes and the quasi-chemical kinetic 
model represents the contribution of the "super-equilibrium" levels of atomic 
oxygen and hydroxyl radical to the total N0 X . Blazowski et al. (ref. 42) 
assembled a model to explain the Lipfert correlation by using a quasi- chemical 
kinetic model to which they added a "prompt NO" contribution. The results shown 
in figure 8 indicate that the so-called prompt NO is due to the overshoot of 
free radical concentrations and the CH + N 2 reactions. 


CONCLUDING REMARKS 

The objective of this study has been to assemble a chemical kinetic 
mechanism that can describe the formation of nitric oxide during the combustion 
of hydrocar bon- based fuels. It was shown that the experimental results could 
be explained in terms of a reaction mechanism in which the CH + N 2 reaction 
process has a significant role in the formation of NO expecially in fuel-rich 
mixtures. The proposed mechanism contains many reactions and rate coefficients 
that have not been experimentally verified. Therefore, the mechanism may not 
be unique and further refinements will be required as more kinetic and mechanis- 
tic data become available. The proposed mechanisms, however, do not appear 
unreasonable and the fact that reasonable agreement was achieved between the 
calculated results and the experimental data for two different hydrocarbons, 
as well as for data obtained from a different reactor (the Lewis Research Center 
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prevaporized, premixed flame- tube- type combustor) , suggests that the essential 
features of the nitric oxide formation processes are contained in the model. 


Furthermore, the similarity between the N0 X levels produced in the propane 
and JP-4 experiments, as well as the ability of the propane combustion and 
nitric oxide formation model to predict the jet-stirred JP-4 results and match 
the Lipfert aircraft gas turbine data, suggest that the behavior of hydrocarbon 
mixtures may not be significantly different from the behavior or pure hydro- 
carbons, such as propane. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
December 1 , 1 980 
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TABLE 1.- METHANE OXIDATION SCHEME 


Reaction 


(1) M + CH 4 -*■ CH 3 + H + M 


1.00 x 10 


43 180 1 


(2) H + CH 4 CH 3 + H 2 

(3) O + CH 4 > CH 3 + OH 

(4) OH + CH 4 -*■ CH 3 + H 2 0 


7.23 x 10 1 


4.10 x 10 14 
3.00 x 10 13 


7 580 1 


7 031 1 


3 020 1 


(5) CH 3 + 0 2 > CH 2 0 + OH 


1.70 x 10' 


7 046 1 


(6) CH 3 + O -*■ CH 2 0 + H 


1.30 x 10 1 


1 006 1 


(7) CH 3 + CH 3 -»■ C 2 H 6 


6.00 x 10 


(8) H + C 2 H 6 > H 2 + C 2 H 5 


1.30 x 10 


4 715 


(9) 0 + C 2 H 6 > OH + C 2 H 5 1.80 x 10 13 

(10) OH + C 2 H 6 H 2 0 + C 2 H 5 6.30 x 10 13 

(11) H + C 2 H 5 * 2CH 3 1.00 x 10 13 


3 070 


1 810 


(12) H + C 2 H 5 + C 2 H 4 + H 2 4.80 x 10 13 

(13) M + C 2 H 5 + C 2 H 4 + H + M 6.80 x 10 17 

(14) H + C 2 H 4 * C 2 H 3 + H 2 1.10 x 10 1 4 


16 004 1 


4 279 


(15) O + C 2 H 4 CH 3 + HCO 

(16) 0 + C 2 H 4 * ch 2 + ch 2 o 


2.50 x 10 


2.26 x IQ 13 


2 516 


1 359 


a The parameters A, n, and C refer to the rate coefficient 
equation k = AT n exp(-C/T) . The rate coefficient units are sec - ^ 
unimolecular reactions, cm 3 /mol-sec for bimolecular reactions, and 
cm*Vmol 3 -sec for termolecular reactions. 


TABLE 1 , 

Continued 



Reaction 

A 

n 

c 


(a) 

(a) 

(a) 

OH + C 2 H 4 + C 2 H 3 + H 2 0 

— i 

• 

o 

o 

X 

o 

0 

1 761 

M + C 2 H 3 C 2 H 2 + H + M 

7.94 x 10 14 

0 

15 853 

H + C 2 H 2 + C 2 H + H 2 

7.76 

3.2 

250 

0 + c 2 h 2 + CH 2 + CO 

5.20 x TO 13 

0 

1 862 

OH + C 2 H 2 * C 2 H + H 2 0 

6.00 x TO 12 

0 

3 523 

C 2 H + 0 2 * HCO + CO 

1.00 x 10 13 

0 

3 523 

C 2 H + O + CO + CH 

o 

X 

_ 1 
o 

CO 

0 

1 585 

CH 2 + 0 2 > HCO + OH 

1.00 X 10 14 

0 

1 862 

CH + 0 2 > HCO + O 

1.00 X 10 13 

0 

0 


(26) CH 3 + H * CH 2 + H 2 

(27) CH 3 + 0 > CH 2 + OH 

(28) CH 3 + OH * CH 2 + H 2 0 

(29) H + CH 2 0 HCO + H 2 

(30) O + CH 2 0 ■+ HCO + OH 

(31) OH + CH 2 0 HCO + H 2 0 

(32) M + HCO + H + CO + M 


1.80 x 10 13 

1.00 x TO 13 

6.00 x 10 10 

3.00 x io 13 
1.20 x 10 14 

6.00 x TO 13 

5.00 x 10 14 


.7 1 510 
.5 8 556 
.7 1 006 
0 2 100 
0 2 200 
0 500 
0 9 562 


a The parameters A, n, and C refer to the rate coefficient 
equation k = AT n exp(-C/T). The rate coefficient units are sec ” 1 
unimolecular reactions, cmVmol-sec for bimolecular reactions, and 
cmVntol 2 -sec for termolecular reactions. 




TABLE 1 Concluded 


Reaction 

A 

n 

c 

Ref. 


(a) 

(a) 

(a) 


(33) H + HCO * H 2 + CO 

1 . 20 x 1 0 1 4 

0 

2 500 

20 

(34) 0 + HCO ^ OH + CO 

1.80 x 10 13 

0 

2 500 

20 

(35) OH + HCO + H 2 0 + CO 

6.00 x 10 12 

0 

0 

20 

(36) OH + 00 + C0 2 + H 

CM 

O 

X 

o 

o 

• 

0 

4 026 

26 

(37) CO + 0 + M * C0 2 + M 

2.50 x 10 15 

0 

2 200 

20 

(38) H + 0 2 OH + 0 

1.22 x 10 17 

-.91 

8 369 

27 

(39) 0 + H 2 OH + H 

2.07 x 10 14 

0 

920 

27 

(40) OH + H 2 ^ H 2 0 + H 

5.20 x 10 13 

0 

3 271 

26 

(41) OH + OH ^ H 2 0 + 0 

5.50 x 10 13 

0 

3 523 

26 

(42) H+H+M>H 2 +M 


1 



M = H 2 0 

5.00 x 10 18 

-1.15 

0 

28 

M = All others 

\ k(M = H 2 0) 



29 

(43) H + OH + M + H 2 0 + M 





M = H 2 0 

io 

** 

o 

X 

o 

0 

-252 

28 

M = All others 

l k(M = H 2 0) 



29 

(44) H + O + M + OH + M 

3.60 x 10 18 

-1.00 

22 

20 


a The parameters A, n, and C refer to the rate coefficient 
equation k = AT n exp(-C/T). The rate coefficient units are sec -1 for 
unimolecular reactions, cm 3 /mol-sec for bimolecular reactions, and 
cm^/mol^sec for termolecular reactions. 
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TABLE 2.- PROPANE REACTION SCHEME 


Reaction 

A 

(a) 

n 

(a) 

C 

(a) 

Ref. 

(45) 

M + C 3 H 8 + CH 3 + C 2 H 5 + M 

5.00 x TO 15 

0 

32 713 

3 

(46) 

CH 3 + C 3 H 8 -*• CH 4 + n-C 3 H 7 

2.00 x 10 13 

0 

5 184 

30 

(47) 

CH 3 + C 3 H 8 +• CH 4 + i-C 3 H 7 

2.00 x 10 13 

0 

5 184 

30 

(48) 

H + C 3 H 8 + H 2 + n-C 3 H 7 

6.30 x 10 13 

0 

4 026 

30 

(49) 

H + C 3 H 8 H 2 + i~C 3 H 7 

6.30 x 10 13 

0 

4 026 

30 

(50) 

0 + C 3 H 8 •> OH + n-C 3 H 7 

5.00 x TO 13 

0 

5 033 

30 

(51) 

0 + C 3 H 8 > OH + i-C 3 H 7 

5.00 x 10 13 

0 

5 033 

30 

(52) 

OH + C 3 H 8 ->• H 2 0 + n-C 3 H 7 

1.60 x 10 14 

0 

1 580 

30 

(53) 

OH + C 3 H 8 + H 2 0 + i-C 3 H 7 

1.60 x lo 14 

0 

1 580 

30 

(54) 

n-C 3 H 7 C 2 H 4 + CH 3 

4.00 x 10 13 

0 

16 658 

30 

(55) 

i-C 3 H 7 -> C 2 H 4 + CH 3 

• 

o 

o 

X 

. — 1 

o 

to 

0 

17 363 

30 

(56) 

n-C 3 H 7 -*• C 3 H 6 + H 

6.30 x 10 13 

0 

19 124 

30 

(57) 

i-C 3 H 7 > C 3 H 6 + H 

2.00 x 10 14 

0 

20 785 

30 

(58) 

0 + c 3 h 6 > c 2 h 4 + CH 2 0 

1.00 x TO 13 

0 

0 

30 

a The 

parameters A, n, and C 

refer to the 

rate 

coeff ici 

ent 


equation k = AT n exp(-C/T) . The rate coefficient units are sec - ^ 
for unimolecular reactions, cm 3 /mol~sec for bimolecular reactions, and 
cm6/rool 2 -sec for termolecular reactions. 
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TABLE 3.- HYDROCARBON FRAGMENT AND NITRIC OXIDE FORMATION REACTION SCHEME 


Reaction 

A 

n 

c 

Ref. 


(a) 

(a) 

(a) 


(59) CH 2 + H -*■ CH + H 2 

2.90 x 10 11 

0.70 

13 085 

25 

(60) CH 2 + 0 * CH + OH 

3.20 x 10 11 

.50 

13 085 

25 

(61) CH 2 + OH + CH + H 2 0 

5.00 x 10 11 

.50 

3 020 

25 

(62) CH + C0 2 > HCO + CO 

1.00 x 10 10 

.50 

3 020 

25 

(63) CH + H > C + H 2 

6.40 x 10 11 

.70 

1 006 

25 

(64) CH + OH * C + H 2 0 

5.00 x 10 11 

.50 

5 033 

25 

(65) C + 0 2 + CO + 0 

5.00 x 10 11 

.50 

2 013 

25 

(66) CH + N 2 HCN + N 

4.00 x 10 11 

0 

6 844 

12 


O 

r— 

O 

X 

o 

o 

00 

0 

6 844 

This study 

(67) CN + H 2 > HCN + H 

6.00 X lo 13 

0 

2 669 

31 

(68) 0 + HCN * CN + OH 

1.40 x 10 11 

.70 

8 505 

32 

(69) OH + HCN * CN + H 2 0 

2.00 x 10 11 

.60 

2 516 

25 

(70) CN + C0 2 + NCO + CO 

3.70 x TO 12 

0 

0 

10 

(71) CN + 0 2 * NCO + 0 

3.20 x 10 13 

0 

505 

31 

(72) H + NCO * NH + CO 

o 

o 

X 

o 

u> 

0 

0 

33 

(73) 0 + NCO + NO + CO 

2.00 x io 13 

0 

0 

33 

(74) N + NCO * N 2 + CO 

1 .00 x io 13 

0 

0 

3 

(75) CH + NO > N + HCO 

• 

o 

o 

X 

o 
— 1 

0 

0 

24 

(76) CH + NO -*■ 0 + HCN 

1.00 x io 13 

0 

0 

24 


a The parameters A, n, and C refer to the rate coefficient 
equation k = AT n exp(C/T) . The rate coefficient units are sec -1 for 
unimolecular reactions , cm*Vmol-sec for bimolecular reactions, and 
cm^/mol 2 -sec for termolecular reactions. 




TABLE 3.- Concluded 


Reaction 

A 

(a) 

n 

(a) 

c 

(a) 

Ref. 

(77) 

0 + N 2 -*• NO + N 

7.60 x 10 13 

0 

38 000 

34 

(78) 

N + 0 2 * NO + 0 

6.40 x 10 9 

1.0 

3 145 

34 

(79) 

H + NO + N + OH 

1.34 x 10 14 

0 

24 760 

35 

(80) 

H + NH + N + H 2 

CM 

O 

X 

O 

O 

.7 

956 

25 

(81) 

0 + NH * N + OH 

6.30 x 10 11 

.5 

4 026 

36 

(82) 

0 + NH NO + H 

6.30 x 10 11 

.5 

0 

36 

(83) 

OH + NH -»• N + H 2 0 

5.00 x 10 11 

.5 

J 006 

25 

(84) 

n 2 o -*■ n 2 + 0 

(b) 



34 

(85) 

n 2 o + co -*• n 2 + co 2 

2.10 x 10 11 

0 

8 757 

37 

(86) 

H + N 2 0 ^ NH + NO 

1.00 x 10 11 

.5 

15 100 

25 

( 87 ) 

H + N 2 0 > N 2 + OH 

7.60 x 10 13 

0 

7 600 

34 

(88) 

o + n 2 o n 2 + o 2 

o 

r— 

X 

o 

o 

0 

14 092 

34 

(89) 

0 + N 2 0 > 2NO 

o 

X 

o 

o 

• 

0 

14 092 

34 


a The parameters A, n, and C refer to the rate coefficient 
equation k = AT n exp(-C/T) . The rate coefficient units are sec*" 1 for 
unimolecular reactions, cm 3 /niol“Sec for bimolecular reactions, and 
cm 6 /mol 2 -sec for termolecular reactions. 

b 1.3 x 10 11 exp(-30 000/T) 

kg 4 = — , where [mJ is the total 

1 + 2.6 x TO -4 LM]" 1 exp (-1 000/T) 
concentration in the system. 
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OF PARAMETRIC STUDIES 


NO x , ppm, as function of fuel-air 
equivalence ratio, 4> 

Methane (300 K) Propane (455 K) 


1 'I 

4 > = 0.8 <t> = 1.0 <t> = 1.2 4 ) = 0.8 4 > = 1.0 4 ) = 1.3 


11 115 I 145 50 205 230 



10 


50 


30 


74 


60 
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oxide formation model to the analysis of NO x emission data reported for various 
aircraft gas turbines showed the contribution of the various nitric oxide forming 
processes to the total NO x formed. 
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